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FOREWORD 


This is a report on Radiation Modeling in the Detection of 
Atmospheric Pollutants. The work was done at the NASA-Langley 
Research Center, during the Summer and Fall of 1972, under the 
Contract No. NAS1-9434-53 . The contract was monitored by 
Dr. Henry G. Reichle. 
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APPROPRIATE LINE PROFILES FOR RADIATION MODELING 


IN THE DETECTION OF ATMOSPHERIC POLLUTANTS 

by Surendra N. Tiwari 
Thermal Engineering Department 
Old Dominion University , Norfolk, Virginia 23508 


SUMMARY 

Absorption by Lorentz, Doppler, and Voight lines are compared 
for a range of atmospheric parameters. It is found that, for the 
intermediate path lengths, the use of the combined Lorentz-Doppler 
(Voight) profile is essential in calculating the atmospheric trans- 
mittance. A brief review of band models, to approximate the absorption 
over certain frequency interval, is presented. Expressions for total 
radiative energy emergent from the atmosphere are given which, with 
appropriate line or band models, can be used to reduce the data obtained 
from radiation measurement by an instrument mounted on an aircraft or 
a satellite. By employing the inversion procedure, the concentration 
of atmospheric pollutants can be obtained from the measured data. 



I . INTRODUCTION 


The study of infrared radiative transmission in actual atmospheres 
requires a detailed knowledge of the atmospheric constituents which absorb 
and emit significantly in the infrared. Specifically this means identi- 
fication of infrared band of major constituents and evaluation of the 
line parameters of these bands. The line parameters depend upon the temp- 
erature, pressure, and concentration of the absorbing molecules and, in 
general, these quantities vary continuously along a nonhomogeneous path 
through the atmosphere. 

With the availability of high resolution spectrometers, it is now 
possible to determine the line positions, intensities, and half widths of 
spectral lines of major atmospheric constituents quite accurately [1-3]. 
Once these fundamental molecular properties are known then it simply 
becomes a matter of employing either a suitable line or band model to 
evaluate the spectral absorption coefficient and subsequently find the 
atmospheric transmittance. 

Because of the complexity of the vibration-rotation spectrum, the 
radiative transmittance is usually calculated by employing the molecular 
band models for homogeneous atmosphere. Through certain scaling approxi- 
mation, such as Curtis-Godson approximation, the solution is then extended 
to the case of nonhomogeneous atmosphere. It should be realized, however, 
that most band models are nothing but interpolation formulas which at 
best correlate the laboratory measured data simulated under homogeneous 
conditions. Consequently these models have limited range of applica- 
bility in pressure, temperature and absorber thickness. While band 



models are very useful in treating radiative transfer problems arising 
in the various areas of gas dynamics [4,5], their validity is restricted 
for actual atmospheric applications [6]. 

The purpose of the present paper is to investigate the feasibility 
of appropriate line profiles for radiation modeling in the detection of 
atmospheric pollutants. Discussions regarding appropriate line profiles 
are presented in Section II. In Section III, integrated absorptions by 
individual lines are calculated by employing the Lorentz, Doppler, and 
combined Lorentz-Doppler (Voight) profiles, and various results are 
compared for a range of atmospheric parameters. A review of different 
band models and band absorptance correlations is presented in Section IV, 
and expressions for upwelling atmospheric radiation are given in Section V. 

The present formulation is based on the assumption of homogeneous 
atmosphere. However, by dividing the atmosphere in appropriate number of 
isothermal layers, the analyses could be extended to the case of non- 
homogeneous atmosphere. 



II. SPECTRAL LINE PROFILE 


In order to describe the infrared absorption characteristics 
of a radiating molecule it is necessary to consider the variation of 
the spectral absorption coefficient for a single line. In general, 

for a single line centered at the wave number to. , this is expressed 
as 


K . 



£ t<“- V • 


( 2 . 1 ) 


where S . 

J 


is the intensity of the j th spectral line and is given by 


S. 
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r 

= J_oo Vj d(O) - w.) . 


( 2 . 2 ) 


The line intensity may be described in terms of the molecular number 
density and Einstein coefficients, i.e., it depends upon the transition 
probabilities between the initial and final states and upon the popul- 
ations of these states. For a perfect gas it may be shown that S 

j 

is a function solely of temperature. The quantity f.(to, y.) f s the 
line shape factor for the jth spectral line. It is a function of the wave 
number to and the line half width y. and is normalized on (to - to.) 

■J j 

such that 
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( 2 . 3 ) 


Several approximate line profiles have been described in the 
literature. Most commonly used profiles are rectangular, triangular, 
Lorentz, Doppler, or Voight (combined Lorentz and Doppler) profiles. 
Lorentz, Doppler, and Voight profiles are of special interest in the 
atmospheric studies. 



The line profile usually employed for studies of infrared 
radiative transfer in the earth's atmosphere is the Lorentz pressure 
broadened line shape for which the shape factor is given by 

f(oj, 0 ) , y ) = y T /{ tt [ (co-co. ) 2 + y t 2 ]} , (2.4) 

J L. Li J Li 

where y l is the Lorentz half-width of the jth spectral line, 

(i.e., half the total line— width at half— intensity) . Variation of 
the spectral absorption coefficient for the Lorentz lines may now be 
expressed as 

Vj = s j " W j )2 + Y L ]} • (2.5) 

From simple kinetic theory, the Lorentz line half-width may 
be shown to vary with pressure and temperature as 

Y l ^ (P//T ) , (2.6) 

where P is the effective pressure and T is the absolute temperature. 
Atmospheric gases which can absorb in the visible and infrared regions 
have relatively small concentrations. Thus, the pressure and the effec- 
tive pressure are not significantly different. 

More detailed quantum mechanical calculations also confirm the 
linear dependency of y l upon pressure, but indicate that the inverse 
square-root variation with temperature is often true only in the wing 
regions (large values of j) of vibration-rotation bands. Considering 
CO 2 example, Yamamoto et al. [7] have shown that the temperature 

dependency of the Lorentz line half-width may be described by y . a, T~ n j , 
and that n^ approaches 0.75 for small j , decreases with increasing -j 
to approximately 0.3, and then increases with a further increase in -j to 
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the kinetic theory value of 0.5. Ely and McCubbin [8] and Tubbs and 
Williams [9] have shown that a more appropriate value for n^ is unity 
for the lines in the immediate vicinity of the band center. 

For Doppler broadened lines, the shape factor is given by the 
expression [10] 

f(w, ay y D ) = (1/Y D ) (in 2 /tt) 1/2 exp [-(oj-gO 2 (£n2/y 2 ) ] , (2.7) 

where y^ is the Doppler half-width and is given by 

Y d - (w /c) (2 k Tin 2/m) 1/2 , (2.8) 

where c is the speed of light, k is the Boltzmann's constant, and 
m is the molecular mass. 

Doppler broadening is associated with the thermal motion of 
molecules. From Eq. (2.8) it is clear that the Doppler width depends 
not only on temperature but also on molecular mass and the location of 
the line center. For certain atmospheric conditions, therefore, the 
Doppler and Lorentz widths may become equally important for a particular 
molecule radiating at a specific frequency. For comparable intensities 
and half-widths, however, the Doppler line has more absorption near the 
center and less in the wings than the Lorentz line (Fig. 1) . 

For radiative transfer analyses involving gases at low pressures 
(upper atmospheric conditions) it becomes imperative to incorporate the 
combined influence of the Lorentz and the Doppler broadening. The shape 
factor for the combined profile is given by [10,11] 

r co 

f(a,v) = (a/TTY D )(£n 2 /tt) 1/2 J (exp(-t 2 ) / [a 2 + (v-t) 2 ]}dt , (2.9) 

^ —00 

where 

3 = ( W (£n 2)1/2 ’ (2 . 10) 
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Fig. 1. Comparison of Lorentz and Doppler line profiles 


having equal line intensities 
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v = [ (to - to )/y ](£n 2) ±/z = (03 - w.)(a/y T ). 

J u J E 


( 2 . 11 ) 


Physically, a represents the ratio of Lorentz to Doppler width and v 
is the wave number scaling factor. 

The appropriate absorption coefficient for the combined profile 
is obtained from Eq. (2.1) and (2.9) as 


/ OO 

(exp(-t 2 )/ [a 2 + (v-t) 2 ]}dt) = < IC(a,v), 

O 


( 2 . 12 ) 


where 


K o = (SjVC^/TO 172 = (S //F)(a/Y L ) . (2.13) 

The function K(a,v) defined in Eq. (2.12) is called the Voight function 
and the line profile expressed by Eq. (2.9) is called the Voight profile. 

Equation (2.12) can be expressed in the alternate form, derived 
initially by Reiche [12], as 

< W j /<0 = K ( a > v ) = (1/ v^TT ) f exp [-at - (t 2 /4) ] cos (vt) dt . (2.14) 

By expanding the integral in (2.14) in a Taylor series, Plass and Fivel 
[13] derived a general expression for K(a,v) in powers of a as 

OO 

K(a,v) = [exp(a 2 -v 2 ) cos(2va)-(2//ir) £ (n ! ) " 1 F (n) (v)a n S in(mr/2) ] , (2.15) 

n=l 

where F^ ^ (v) is the nth derivative of the function 

r 00 r 

F(v) = / exp (-t 2 ) sin(2vt) dt = exp(-v 2 ) j exp(t 2 ) dt . (2.16) 


For I w I » Y l and | co - oi | » y d Eq . (2.15) reduces to 


asymptotic form as 


K(a, v) = (a/v 2 /rF) [1 + (| - a 2 ) + (~ - 5a 2 + a 4 ) + ] . (2.17) 



For (v/a) » 1 , Born [14] has given an asymptotic expansion for 
K(a,v) as 

K(a,v) = {a/ [ (a 2 + v 2 ) SH ]}[1 + (3/2 v 2 )] . (2.18) 

For small values of the parameter a , Eq. (2.15) may be expressed 
by [15] 

K(a,v) = exp(-v 2 )-(2a//rr) [l-2vF(v) ] +a 2 (l-2v 2 ) exp(--v 2 ) 

- (2a 3 //F)[ |(l-v 2 ) - 2v(l - | v 2 )F(v) ] 

+ a 4 ("I - 2v 2 + ~ v 4 ) exp(-v 2 ) . (2.19) 

By combining Eqs . (2.1) and (2.7), the absorption coefficient 
for the Doppler broadened lines may now be expressed as 

Vj = K q exp (-v 2 ) . (2.20) 

It should be noted that as the pressure becomes small, y 

L 

approaches zero and Eq. (2.12) reduces to Eq. (2.20). On the other 
hand, for large pressures, the quantity (y /y ) becomes large and 
Eq. (2.12) reduces to the Lorentzian case, Eq. (2.5). In other words, 
the Voight profile assumes the Lorentzian shape in the limit of large 
v and reduces to the Doppler profile for small a . 

Extensive reviews of the literature on the Voight function and 
some accurate approximations to the Voight profile are given in 
[10, 11, 16-18]. Armstrong [17] has presented a review of the mathe- 
matical properties of the Voight function and has discussed the various 
computational methods which make efficient and economical use of computer 
time and storage. With the aid of computer programs developed by 
Armstrong [17], Young [18, 19] or Hummer [20, 21], it is now possible to 



calculate the Voight function, for an extended range of parameters, 
to an accuracy of better than six significant figures. 

A simple closed form approximation to the Voight profile 
that is valid over a useful range of parameters is given, in terms 
of the present nomenclature, by [22-24] 

K wj /K 0)o = [1 ~ ( Y L /Y V )J e xp{- H.088 [(w-a)J/Y v ] 2 } 


+ (Y L /Y V ) /{1 + 16 [(oho )/ Y 12} , 


v- 

where Voight half-width is expressed in terms of Y and Y as 

L 'D 


Y v = (Y L /2) + [( Y 2/4) + ] 1/2 


(2 


and 


(2 


K o)o S j/^Y V [1.065 + 0.447(y l /y v ) + 0.058( Y / Yv ) 2 ] > 

L v 

This form is very convenient for numerical computation and it matches 

Che Voight profile within 5 per cent under worst conditions. General] 

the error is within 3 per cent, with maximum errors occurring near 
zero pressures. 


• 21 ) 

22 ) 

23) 



III. ABSORPTION BY INDIVIDUAL LINES 


The radiative transmittance at a single wave number may, in 
general, be expressed as 


\ - ex p ( - / 

* / o 


X 


S f dX) , 


(3.1) 


where X is the mass of the absorbing gas per unit area and is equal 
to the density times the thickness of the absorbing layer. For a 
nonhomogeneous atmosphere, X is given by the relation 


X = 


£ 


p a 
a 


(3.2) 


where £ is the length measured along the direction of the path which 
makes angle 0 with the vertical, and p is the density of the 

3 . 

absorbing gas. For a homogeneous path, Eq. (3.1) reduces to 


Vj = ex P( -S j f x > = exp(-K X) . (3.3) 

For n independent lines within the wave number range of interest, 

the absorption coefficient is given by 

n 


< = V K . , 

U) W] 

3=1 


(3.4) 


and consequently the transmittance is expressed as 

n 


T 


= n 

3=1 


T . . 


(3.5) 


The total (integrated) absorption for a single line is given by 


/ CO 

(1-T .) d (to - w. ) 
3 


(3.6) 


where again ul represents the wave number at the line center of the 
jth spectral line. In astronomical and meteorological literature, the 
relation given by Eq. (3.6) is usually referred to as the equivalent 



width of a single line. 


1. Single Line with Lorentz Shape 

For homogeneous atmospheric path, the absorption of an isolated 
spectral line of Lorentz shape, A^, may be obtained by combining 
Eqs. (2.5), (3.3) and (3.6) as 


f°° 

K = I { 1-exp (-S . y X/tt[ (u)-ti). ) 2 + y 5 ] ) } d(oj-o).) . 

L J-oo 3 L 1 L J 


(3.7) 


In this equation, the extension of the limits to infinity is physically 
realistic since the absorption beyond a certain wave number from the 
line center is zero. By defining nondimensional quantities 


x = S.X/2-ny , 
J L 


y = (w-oO/y , 
3 ^ 


Eq. (3.7) can be written as 

\ = 2Y L [ ^ “ exp[-2x/(y z + 1) ] }dy • 

J n 


(3.8a) 


(3.8b) 


(3.9) 


An exact solution of this equation is found in terms of the Landenberg- 
Reiche function, L(x) , as [25, 26] 


\ = 27TY L L< - X ^ = 27 Fy l {x exp(-x)[I Q (x) + I 1 (x)]} , 


(3.10) 


where I q and 1^ are the Bessel functions of imaginary arguments. 

Some values of the function L(x) are tabulated in references [10,27]. 
By multiplying Eq. (2.5) by the mass of the absorbing gas, X, it can be 
shown that the quantity x represents one-half the optical path at the 
line center. 

Equation (3.9) possesses two well known asymptotic forms. For 
small x, expanding the exponential in Eq. (3.9) in series and retaining 



only the first two terms, results in 


\ = 4 Y l x j (y 2 + 1) 1 dy = 2 tty l x = S^X . (3.11) 

This is called the linear limit or weak line approximation. A line 
is called weak when K^X << 1 for all frequencies within the line. 

In this limit, absorption is independent of the pressure and is a 
linear function of the amount of the absorbing gas and the line 
intensity. According to Plass [28,29], this approximation is accurate 
within p per cent if x < 0.02p. 

For large values of x, Eq. (3.9) reduces to 


Al-^l 


[1 - exp (-2x/y 2 ) ] dy , 


(3.12) 


which may be expressed in an alternate form as 
Aj^ = 2y L (2x) 1//2 / {[1 - exp(-z 2 ) ] / z 2 }dz , 

where z 2 = 2x/y 2 . Integration of Eq. (3.13) yields 


(3.13) 


A l = 2y l (2ttx) 1/2 = 2(S Y l X) 1/2 


(3.14) 


This is the square-root limit or the strong line approximation and is 
accurate within p per cent if x _> (12.5/p). Since the line width 
varies linearly with the pressure, absorption in this limit varies as 
the square root of both the pressure and the amount of the absorbing 
medium. Almost all incident radiation, in this limit, is absorbed at 
frequencies within the total half width of the line center. 

Equations (3.11) and (3.14) can also be obtained from Eq. (3.10) 
by expanding the modified Bessel functions for small and large values of 
the argument respectively [26,27]. For an extensive physical explana- 
tion of the weak and strong line approximations, see Penner [10] and 
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Goody [27]. The solution expressed by Eq. (3.10) is illustrated in 
Figure 2 along with the limiting solutions. 

Three approximate relations for Eq. (3—9) are suggested in 


references [30,31] as 


\ = (2tty l ) (2 x/tt) 1/2 {1 - exp[-(TTx/2) 1/2 ]} , 

(3.15a) 

A l = (2tty l ) (2 x/tt) 1//2 [1 - exp(-7Tx/2) ] X ^ 2 , 

(3.15b) 

A l = ( 2 tty l )/[1 + (ttx/2)] 1/2 . 

(3.15c) 


Equations (3.15) reduce to the correct asymptotic limits, and their 
solutions agree well with the exact solution, Eq. (3.10), over the 

entire range of x . Results of Eq . (3.15c) are also illustrated in 
Figure 2. 

2. Single Line with Doppler Profile 


Doppler broadening becomes important under upper atmospheric 
conditions. The absorption of a Doppler broadened line, A^, is 
obtained by combining Eqs. (2.20), (3.3), and (3.6) as 

/■GO 

Ajj = [Y D /(£n2) 1 2 ] I (1 - exp[-x exp(-v 2 )]}dv , (3.16) 

—OO U 


where 


X D = K o X = f(S j /Y D )(£n2/TT) 1/2 ] X 


represents the optical path at the line center. By expanding the 


(3.17) 


exponential in an infinite series, Eq. (3.16) can be written as [10] 
Ap = (fT/2.n2) 1/2 Y;d f(x D ) , 


(3.18a) 



12 *y L 



Fig. 2. Absorption by a single line of Lorentz shape 




16 


where 00 

f(x D ) = H ("D n x D n+1 /[(n + 1)! (n + 1) 1/2 ] . (3.18b) 

n=o 

This may be expressed in the nondimensional form 

Ajj = A D /[2y D /(£n2) 1/2 ] = (/fr/2) f(x D > . (3.19) 

For 0.1 £ x^ <_ 30 , Eq. (3.19) may quite accurately be 
approximated by 

Ajj 2L C^/ 2 ) x d exp[-0.5 (x D ) 1/2 ] . (3.20) 

As x D - o (linear limit), Eq . (3.19) reduces to 

A d = (v^F/2) x q , (3.21a) 

or 

Ap = S j X ’ (3.21b) 

This is identical to the result for the Lorentz profile in the linear 
limit (Eq. 3.11). Moreover, absorption by any line shape must reduce 
to this expression in the linear limit. This is because the linear 
limit can be obtained, independent of the line shape, directly from 
Eq. (3.6). 

For large values of x D , Eq. (3.19) may be expressed by an 
asymptotic expansion of the form 

I/O 00 

S = X D ) (l-[l/2r(l/2)] £ r (n) (1) r (n-f)/[n!(£n x D ) n ]} , (3.22) 

(ft) 

where F (1) is the nth derivative of the gamma function evaluated 
at unit argument. In the limit of large x^ , this is further approxi- 
mated by 

1/2 


% - (£n V 


(3.23) 



It is important to note that in the limit of large x D , the 
absorption increases very slowly with the amount of absorbing gas. 

This is because in the large path length limit the central portion 
of the line becomes opaque and absorption occurs solely in the line 
wings. Since the Doppler line shape drops off exponentially in the 
line wings, it is not possible to absorb much additional radiation 
in the wings by increasing the number of molecules. 

The solutions expressed by Eqs. (3.19) and (3.20), as well as 
the limiting solutions, are illustrated in Fig. 3. 

3* Lines with Combined Doppler and Lorentz Profile 

For conditions where both the Doppler and the Lorentz broaden- 
ings are important, the total absorption, , is obtained by combining 
Eqs. (2.12) and (3.6) as 

/ OO 

U - exp [-x K(a,v) ] }dv . (3.24) 

— OO 

It can easily be shown that Eq. (3.24) reduces to the Lorentzian case 
for large values of a and to the Doppler case in the limit of small a. 

Numerical solution of Eq. (3.24) may be obtained by employing 
the appropriate expression for the Voight function, K(a,v). Using the 
numerical procedure described by Young [18,19], Eq. (3.24) was solved 
for a range of the parameter a, and the results are illustrated in Fig. 4 
Similar results were also obtained by Jansson and Korb [32], who employed 
a somewhat modified version of Armstrong's [17] Fortran IV program. These 
results are in general agreement with the results of Penner and Kavanagh 
[16]. As would be expected, for a = o, the results correspond to the 
case of a pure Dopper profile, while for a > 10 they correspond to the 



Strong Line 



Approximate, Eq. (3.20) 






Lorentzian shape. 

If the approximate relation for a Voight Profile, as given by 
Eq. (2.21), is employed, then the total line absorption, Ay , is 
expressed by 

r oo 

Ay = Ay/[2y D /(£n2) 1 2 ] = j I {l - exp[-x D W(a,v)]}dv , (3.25) 

—OO 

where 

W(a,v) / [ (£n2/ /rT)F(a) ] = [ 1— (Y L /Y y ) ] exp{-11.088[(v/a) ( y l /y v > ] 2 } 

+ (Y l /Y v )/{1 + 16[(v/a)(Y L /Y v )] 2 } , 

F(a) = (a[1.065/(Y L /Y v ) + 0.447 + 0.058 (y^)])" 1 , 

(Y L /Y V ) = {0.5 + [0.25 + (£n2/a 2 )] 1/2 } _1 . 

The numerical solution of Eq. (3.25) is straightforward. In Fig. 5, 
some results of Eq. (3.25) are compared with the solutions obtained by 
using the actual Voight profile. Except for very small values of the 
parameter a , the agreement between the two solutions is quite satis- 
factory. 

4. Comparison of Integrated Line Absorption 

The errors, in the integrated line absorption, encountered in 
using the Lorentz shape to approximate the Voight profile, are illus- 
trated in Fig. 6 for a range of the parameter a . x As would be expected, 
maximum errors occur for intermediate path lengths and for lower values 
of a . The errors, however, are not significant in the two limiting 
cases of small and large path lengths. As discussed earlier, for small 
path lengths (linear limit) , the line absorption is independent of the 
line shapes. In the limit of large path lengths, however, the central 




a = 10 





x 100 
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portion of the line becomes opaque and absorption occurs only in the 
wing regions. Since the Voight profile is essentially Lorentzian in 
the line wings, the error encountered in the large path length limit 
becomes insignificant. 

The errors resulting from using the Doppler line profile to 
approximate the absorption by the Voight line profile are illustrated 
in Fig. 7. Since a = o corresponds to the case of pure Doppler 
absorption, the errors are expected to be higher for larger values of a . 
Maximum errors, in this case, are found in the large path length limit. 

This is because, in this limit, the absorption occurs essentially in the 
line wings and the Voight profile in the wing regions is Lorentzian 
rather than of Doppler shape. 

It should be emphasized that for cases of intermediate path lengths 
and for moderate values of a (0.1 < a < 10) the use of the Voight pro- 
file becomes almost essential. This situation corresponds to the radiative 
transmittance in the earth's troposphere and lower stratosphere. Conse- 
quently, for radiative modeling of the lower atmosphere, consideration must 
be given to the application of combined line profiles. 

In Figs. 8-12, comparisons of results for all three line profiles 
(Lorentz, Doppler, and Voight) are made for a = 0.001, 0.01, 1.0, and 10. 
These figures clearly illustrate the range of validity of the absorption 

intersect at some point, say x^, the location of which increases with 
decreasing a . For a > 2 , and curves do not intersect, and 

Al and curves become identical for almost all path lengths. 


by the three line profiles. In Figs. 8 — 10, the curves for 



)/ArJ x 100 
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Yamada [33] has reported that A^ and A^ curves became identical 
for a > 0.2 which is not in agreement with the present investigation. 

It should be realized, however, that the parameter a is a measure of 
the relative importance of Lorentz versus Doppler effect. For a value 
of a = 0.5 , for example, the ratio of Lorentz to Doppler width, 

(Y l /Y d ) = a//£ n~2 n* 0.6 . In this case, therefore, the Doppler effect 
will be significant especially for the intermediate path lengths and 
integrated absorption could not be expressed solely in terms of the 
absorption by Lorentzian profile. It is quite possible that the way 
Yamada' s results were illustrated, the difference between A^ and A^ 
curves was not apparent in the intermediate path lengths. For small as 
well as large path lengths, A^ and A^ are, of course, identical. 



Log ln ( A) 
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Fig. 9. Comparison of line absorption for a = 0.01 
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IV. REVIEW OF BAND MODELS AND BAND ABSORPTANCE CORRELATIONS 


The absorptions by atmospheric gases, under certain atmospheric 
conditions, result in overlapping spectral lines. Thus, the total 
absorption within certain frequency interval cannot accurately be 
represented simply by summing the absorption by individual lines. 

This is because the absorption in a region of overlapped lines is 
always less than the absorption calculated by considering the contri- 
butions of individual non-overlapping lines. 

The total absorption of a band of overlapping lines strongly 
depends upon the line intensity, the line half-width, and the spacing 
between the lines. In a particular band, the absorption coefficient 
varies very rapidly with the frequency and, therefore, it becomes a 
very difficult and time-consuming task to evaluate the total band 
absorptance by numerical integration over the actual band contour. Con- 
sequently, several approximate band models have been proposed which 
represent absorption from an actual band with reasonable accuracy. 


1. Elsasser Band Model 


The simplest band model that accounts for the line structure is 
the Elsasser model [26] , for which equally spaced lines of equal intensity 
and equal half width are assumed. The expression for the total band 
absorption by Elsasser band model has been developed only for lines 
having the Lorentz line profile. For homogeneous atmosphere, the total 
absorptance of the Elsasser band is found to be [5,27] 


■ (A q /2) 


f 


{l-exp [- 


u sinh(ir3/2) 


cosh(fr$/2) - cos(ttz/2) 


] } dz 


(4.1) 
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where 

A = nd , u = Spy/A , 6 = 4y, /d , z = 4(to-w )/d . 

O O Li O 

In this equation u is the dimensionless pressure path length, 3 is 
the line structure parameter, A q is the band width parameter, n is 
the number of lines in the band, and d is the spacing between the 
lines. Equation (4.1) has a form that is characteristic of all band 
models, namely that the total absorptance may be expressed as 

A = A q A(u,0) , (4.2) 

where A(u,8) is a dimensionless function. Equation (4.1) reduces to 
several limiting forms which are discussed in [5] . By employing the 
scaling approximations mentioned in references [34-37], Eq. (4.1) may 
be used to evaluate the total band absorptance for the case of hon- 
homogeneous atmosphere. 


2 . Statistical (Mayer - Goody) Band Model 


The statistical band model is based upon the assumption that, 
in a given wave number interval, the spectral lines are randomly 
spaced and that the intensity of these lines can be specified by some 
distribution function. For this model, the absorptance over a wave 
, i number interval , D = nd , is given by [27-29] 

A = 1 - (1 - A jjD ) n (4.3) 

where n is the number of lines in interval D , d is the mean 
spacing between the lines, and 


A J; d(S 0 , x ,P) 


-f 

•' n 


A (S,X,p) P (S , S )dS. 
J o 


(4.4) 
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In Eq. (4.4), A^,D is the absorptance of a single line over the entire 
wave number interval D , P(S,S o ) is the normalized probability of 
finding a spectral line with an intensity between S and S + dS, and 
S q is the parametric line intensity that occurs in the intensity dis- 
tribution function. 

By employing one definition of the exponential, it may be shown 
that for large number of lines in the interval D, Eq. (4.3) becomes [29] 


A = 1 - exp (-n A ) , 
1 » u 


n » 10 


(4.5) 


If it is assumed that all lines, within the wave number interval 
D , are equally intense then 

P(S) = 6(S - S ) (4.6) 

o 

and it follows that 

h,D ( V X ' p) ’ A l.D (f V X>P) ' (4 ' 7) 

Now, the expression for A^ Q can be obtained from Section-Ill for 
a particular line shape simply by replacing S by S q . 

The development of the expressions for absorption by a single 
isolated line, in Section-Ill, assumes an infinite wave number interval. 
Consequently, these expressions can be used in Eqs. (4.4) and (4.5) only 
if there is no appreciable absorption by the individual line beyond wave 
number interval D. If this is not the case then the expression for 
A_. D should be obtained by integrating the expression for the line 
absorption over the actual frequency interval D. 

If an exponential distribution of line intensities is assumed, then 


P(S) = exp (-S/5 )/S 
o o 


(4.8) 
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and, for a particular line shape, the line absorption can be calculated 
by using this line intensity distribution. For spectral lines with 
Lorentz shape, the absorption by the statistical band model, in this case, 
is found to be [29] 

A = l-{l-(TT0x o /2)/n(l + 2x o ) 1/2 } n (4.9) 

or 

A = l-exp{-(TT0x o /2)/(l + 2x q ) 1/2 } for n » 10 , (4.10) 

where 

X = S X/27TY 
o o ' L 

Since the lines are assumed to be randomly distributed in the 
statistical model, the absorption by this model is always less than that 
by the Elsasser band model. The advantage of the statistical model is 
that it can be applied to any line shape. 

3. Random Elsasser Band Model 

In actual vibration-rotation bands , lines are arranged neither 
completely at random nor at regular intervals. Within a band, there may 
be a number of strong lines in certain narrow spectral region while in 
other regions only very weak lines may be present. There may also occur 
a superposition of equally intense lines in another spectral region 
of the band. In cases like this, a more accurate representation of band 
absorption is provided by the random Elsasser model, which assumes the 
random superposition of several different Elsasser bands. Each of the 
superposed bands may have different line intensities, half-widths, and 
spacing. As many different Elsasser bands as necessary may be superimposed 
in this model. As the number of superposed Elsasser bands becomes large, 
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the absorption by the entire band approaches that given by the statistical 
band model. For N randomly superposed Elsasser bands, the absorption 
is given by the relation [27,29,38] 

N 

A = i - n [1 - A ,(x.,B.) (4.11) 

i=l 1 1 


where 



is the absorptance of an Elsasser band. 


4. Quasi-Random Band Model 

Quasirandom model is probably the best model to represent the 
absorption of a vibration-rotation band quite accurately. It assumes 
neither a regular nor a random spacing of the spectral lines. The 
essential feature of this model is to divide the wider frequency 
interval of the actual band into much narrower subintervals. In each 
of these subintervals, the spectral lines are assumed to have random 
spacing. In this manner, the model accounts for the actual intensity 
distribution of strong as well as weak spectral lines. The absorption 
of the narrow spectral interval is calculated from the relation of a 
single-line absorption over a finite interval. The absorption of each 
of the n lines in the narrow interval is calculated separately and 
the results are combined by assuming a random position for the lines 
within the interval. The total band absorptance is calculated by averaging 
the results from the smaller intervals. 

The total band absorptance of a quasirandom band model is given 
by the expression [39] 

1 6 

^ = J H A k 

K=1 


(4.12) 
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where 


M 

a = i - n [i 

K i=i 



(4.13) 


In these equations, 6 is the width of the narrow frequency intervals 

and A is the corresponding absorptance of each subinterval. A. ^ 

K 3 >D 

is the absorptance of a single line over the finite interval D and 
M is the number of lines in the frequency interval k . 

The four band models reviewed above are called the narrow band 
models. There are three limits in which approximate expressions for the 
total band absorptance, by these bands, can be obtained. These are 
weak-line approximation, strong-line approximation, and nonoverlapping- 
line approximation. These approximations are discussed in detail in 
reference [40]. 

The narrow band models may be employed to compute the transmittance 
from atmospheric gases under most atmospheric conditions. Before the use 
of these band models can be justified, however, it would be advisable to 
compare them with the line-by-line model especially for gases whose line 
parameters are completely known. 

5 . Wide Band Models 

Aside from the narrow band models discussed above, there are also 
available in literature the so-called wide band models which provide 
correlations that are valid over the entire band pass. Besides possessing 
the conventional linear and square-root limits, these models also possess 
another asymptotic limit which is called the logarithmic limit. Even 
though the use of these models may be restricted for atmospheric appli- 
cations, they do provide a quick and accurate information regarding 
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transmittance of gases at moderately high temperatures. As such, 
these models are quite useful in many engineering applications. 

a) The Box Model : The simplest of the wide band models is the 

box model which was first introduced by Penner [10]. For this model, 

it is assumed that the absorption coefficient, < , is constant over 

CO 

an effective band width Acj . The expression for the total absorptance 
by this model is given as 

A = / [1 - exp (-K X) ] dw = (Acj) [1 - exp(-icX)] 

•'Aw ■ e 

where (Acj)^ is the effective band width, and K is the mean absorption 
coefficient for the interval (Aoo) e • Further discussion and application 
of the box model is available in references [10,31,41] . 

b) The Exponential Wide Band Model : Edwards and his co-workers [4,42] 

have considered various wide band models (Rigid rotator, non-rigid 
rotator, and Arbitrary models) and have concluded that three parameters 
(the mean line intensity to spacing ratio, the mean line-width to spacing 
ratio, and the effective broadening pressure) are necessary for a complete 
description of the band absorption. For complete discussions of these 
models and their limiting forms, one should refer to [4,5,31,42]. By 
using the Curtis-Godson approximation, the exponential wide band model 
can be applied to calculate the band absorption for nonhomogeneous gases. 

6. Band Absorptance Correlations 

Three continuous relations, for the total band absorptance, 
which are valid over all values of path lengths and line structure para- 
meters, are available in literature [5,31,43]. These are. 



(a) Tien and Lowder: 


A = A/A = £n(uf(3){(u + 2)/[u + 2f(3)]} + 1) , 
o 

where 

f (3) = 2.94[1 - exp (-2. 603) ] , 

(b) Goody and Belton: 

A = 2 An{l + [uvf/(u + 43) 1/2 ]} , 

(c) Cess and Tiwari: 

A = 2 £n{l + u/[2 + Ai(l + T7F) ] ) • 

Various limitations as well as applications of these correlations are 
discussed in [5]. 

By employing either an appropriate line or band model, the total 
(integrated) band absorptance can be obtained for a range of atmospheric 
parameters. It would be interesting to compare the band absorptance 
results obtained by using different line and band models. 
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V. EXPRESSIONS FOR UPWELLING ATMOSPHERIC RADIATION 


As shown in Fig. 13, the radiation emergent from the atmosphere. 


E(to) , may be given by the expression [44] 


E (oj) = E g (oj) + E r (w) + E^Coo) + E R(j) (a)) 


(5.1) 


where 

E (w) = thermal radiation emitted by underlying surface and 
Cj 

atmosphere 

E ( to) = incident solar radiation reflected by the surface 
E^ (to) = radiation scattered by a single or multiple scattering 
processes in the atmosphere without having been 
reflected from the surface 

Er^U)) = scattered energy which has undergone a reflection from 
the surface. 

In general, these quantities are functions of surface temperature, 
atmospheric temperature, surface emittance, surface reflectance, sun 
zenith angle, scattering characteristics of particles, and transmittance 
of the atmosphere. 

Neglecting the scattering and solar radiation, the equation of 
radiative transfer for thermal radiation emerging from a plane-parallel 
atmosphere, with underlying surface of emittance e(to) , can be written as 


E(to) = E (oj) = e(w) B(w,T ) t(w,o) 
Ct s 


+ 1 B (a),T(z) ) [dT (co,z)/dz]dz 


(5.2) 


where B(co,T) is the Planck's blackbody function, T g is the surface 
temperature, T(z) is the temperature at altitude z, and t(w,z) is the 
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monochromatic transmittance of the atmosphere (at wave number to ) 
between altitude z and the top of atmosphere. 

The contribution from sunlight reflected from the surface 
becomes significant at shorter wavelengths. This contribution is 
given by the component E R (to) as 

E R (to) = (1 /tt) [1 - e (to) ] cos 0 H g (to)[T(to)] ? (5.3) 

where 0 is the sun zenith angle and 5=1+ f ( 0) . Function 
f(0) = sec 0 for o < 0 < 60° and equals to Ch 0 for 0 > 60° with 
Ch 0 denoting the Chapman function. H g (to) is the sun irradiance on 
top of the atmosphere, and x(to) is the transmission vertically through 
the atmosphere. 

The expression for the transmittance x(to,z) is obtained from 
an appropriate molecular (band or line) model considered in the previous 
sections. In general, the transmittance may be expressed as 

x(to,z) = n [T.(to,z)] = exp [- I £ 6 i (w,z") n ± (z') dz'] (5.4) 

i J o i 

where B^(to,z) is the extinction coefficient for species i which is 
equal to the sum of absorption and scattering coefficients, and n^ is 
the number density of the contributing molecules of species i in the 
optical path. 

A combination of Eqs. (5.1) through (5.4) yields a relation for 
thermal radiation emerging from a plane-parallel atmosphere. This 
equation, with appropriate spectroscopic information, can be used to 
reduce the measured atmospheric data by an inversion procedure. The 
equation may also be used to obtain the concentration of atmospheric 



pollutants, from radiation measurements, provided other governing 
parameters are known. 

In radiation modeling for pollution measurement, the upwelling 
radiation from the nonhomogenous atmosphere is calculated by dividing 
the atmosphere into an appropriate number of sublayers. Each of these 
sublayers is assumed to be homogeneous in species concentration, 
temperature and pressure. 

It should be noted that the derivation of Eq. (5.2) assumes the 
existence of local thermodynamic equilibrium. This assumption, however, 
is not justified for gases like carbon monoxide for the range of atmo- 
spheric pressures and temperatures. Since the nonequilibrium effect 
comes only through the source function [27,45], for gases like CO, the 
Planck function in Eq. (5.2) should be replaced by the nonequilibrium 
source function. The use of Eq. (5.2) is, however, justified for gases 
like CO 2 , Chi,, and H 2 O. 



VI . CONCLUDING REMARKS 


The purpose of this study was to compare the total (integrated) 
absorptions of the three line profiles, Lorentz, Doppler, and Voight, 
for a range of the governing parameters. It is concluded that the Voight 
line profile should be employed in calculating the transmittance for the 
troposphere and lower stratosphere. Line-by-line models are preferred 
over band models for atmospheric applications. The line-by-line model 
is especially convenient in cases where there is interference of other 
species within the spectral region of interest. 

The equation for the radiative transmittance, based on the line- 
by-line model, is used in the expressions for upwelling radiation which 
is received by an aircraft or satellite mounted instrument. From these 
radiation measurements, amounts of atmospheric pollutants are determined 
by employing an appropriate data reduction scheme. The essential para- 
meters, appearing in actual calculation of the amount of each pollutant, 
are: ground temperature, ground emissivity, sun zenith angle, atmospheric 

temperature profile, water vapor distribution, and distribution of inter- 
fering species. Furthermore, within the spectral region of each pollutant, 
there may appear several thousand spectral lines.. Thus, time required for 
data reduction becomes quite large. One way to reduce the data reduction 
time will be to employ an appropriate narrow band model in place of the 
line-by-line model. If the spectral region of interest could be divided 
into a number of convenient narrow spectral intervals, then it is quite 
possible that either the statistical or quasi-random band model would pro- 
vide an accurate description of the entire spectral transmittance. However, 
before these band models can be employed for atmosphere applications, their 



validity must be established by comparing the transmittance of these models 
with the results obtained by using the line-by-line model. 
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LIST OF SYMBOLS 


A 

A 

V 

% 

A E,i 

A. 

J 

A. „ 

3 ,D 

*L 

*L 

A 

o 

\ 

\ 

B(co,T) 

c 

E 

e g 

e r 


J R<j) 


parameter representing the ratio of Lorentz to 

Doppler width, Eq. (2.10) 

total band absorptance, cm 1 

dimensionless band absorptance 

total absorptance of a single Doppler line 

dimensionless A^ 

absorption of a narrow Elsasser band 

total absorption of a single line (jth line) 

absorption of a single line over wave number interval D 

total absorption of a single Lorentz line 

dimensionless A^ 

band width parameter, cm 1 

total absorption of a single Voight line 

dimensionless A^ 

total absorption of a single line given by Eq. (3.25) 
dimensionless A^ 

Planck's function 

speed of light 

upwelling radiative energy 

thermal radiation of ground and atmosphere 

solar radiation reflected from ground 

solar radiation scattered by the atmosphere 

without having been reflected by the ground 

solar radiation reflected from ground outside the field of 

view and scattered by the atmosphere into the field of view 



h 

H 

s 

K 

K(a,v) 

m 

P,P 

s j 

T 

T 

s 

T(z) 


v 


D 


D 


Y t 


'V 


e 

0 

K 


0 ) 


K . 


altitude of an aircraft or a satellite 

sun irradiance at the top of the atmosphere 

Boltzmann constant 

Voight function, Eq. (2.12) 

molecular mass of the absorbing medium 

partial pressure of the absorbing medium 

line intensity or line strength, cm 2 

kinetic temperature, °K 

surface or ground temperature 

temperature of the atmosphere at altitude z 

dimensionless coordinate, u = Spy/A 

wave number scaling factor, Eq. (2.11) 

mass of absorbing gas per unit area 

optical path at the line center, Eq. (3.8a) 

optical path at the line center, Eq. (3.18a) 

optical path at the line center, Eq. (4.10) 

line structure parameter 

line half-width of a spectral line 

Doppler line half -width 

Lorentz line half-width 

Voight line half-width 

surface emittance 

sun zenith angle 

absorption coefficient correlation, Eq. (2.13), 
equilibrium spectral absorption coefficient, cm 
absorption coefficient of jth spectral line 


cm 

l 
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< 

too 

T 


T 

0) 

0) 


to. 

J 


spectral absorption coefficient correlation, Eq. (2.23). 

radiative transmittance 

spectral transmittance 

wave number, cm 1 

wave number at the line center 
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